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Abstract 



The decays of D mesons to Klv and K*lv final states exhibit significant deviations from 
the predictions of heavy-quark symmetry, as one might expect since the strange quark's 
mass is of the same order as the QCD scale. Nonetheless, in order to understand where the 
most significant effects might lie for heavier systems (such as B — > Dlv and B — > D*lv), 
the pattern of these deviations is analyzed from the standpoint of perturbative QCD and 
0(l/m s ) corrections. Two main effects are noted. First, the perturbative QCD corrections 
lead to an overall decrease of predicted rates, which can be understood in terms of production 
of excited kaonic states. Second, 0(l/m s ) effects tend to cancel the perturbative QCD 
corrections in the case of Klv decay, while they have minimal effect in K*lv decay. 
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1 Introduction 

In atomic physics, one can often separate the nuclear and electronic effects from one another. 
Similar progress has been made for strongly interacting particles within the past few years 
|I],@,[3|,f§|- Particles containing a single heavy quark can be thought of as analogues of atoms 
for the strong interactions. The single heavy quark behaves as the nucleus, while the light 
quarks and gluons behave as the electron cloud. To a large extent the degrees of freedom 
of the heavy quark and those of the remaining matter can be discussed separately. The 
resulting description validates some (but not all) older quark-model results in a systematic 
way, and has come to be known as heavy quark symmetry. It represents a limit of quantum 
chromodynamics (QCD) for infinite mass tuq of the heavy quark Q. 

The semileptonic decays of mesons containing a single heavy quark are an ideal labora- 
tory for testing heavy quark symmetry. Both the leading-order results and the 0(1 /tuq) 
corrections to them have been fully worked out. In practice the symmetry is applicable 
only to the semileptonic decays B — > Dlu, B — > D*lu, and possibly to decays involving ex- 
cited D mesons in the final state. Only the b quark (in the B meson) and the charmed quark 
(in the D and D*) are heavy enough compared to the QCD scale of several hundred MeV 
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that one can begin to think of applying heavy quark symmetry. It is possible to describe 
six form factors (two for B — > Dlv and four for B — > D*lv) in terms of a single universal 
function of the variable w = v ■ v', where v and i>' are the four- velocities of the initial 
and final states. We shall refer to £ as the Isgur-Wise function. 

One would like very much to know the size of corrections to the universal Isgur-Wise 
behavior in B semileptonic decays. One benefit of this information would be the improved 
ability to not only measure the b — > c weak coupling accurately J7LH], but to be able to reliably 
estimate the theoretical error. However, up to now it has not been possible to estimate the 
actual magnitude of the corrections. 

Early estimates of semileptonic decays were not restricted to cases in which heavy-quark 
symmetry was guaranteed to work. In particular, many attempts were made to describe the 



semileptonic decays of D mesons within the framework of various quark models P,|10|,|iT 
|12|) |13|, |14l . The main shortcoming of these descriptions was that they predicted rates for 
D — > K*lv which were too large. An attempt was made |TS| to isolate the source of the 
discrepancy by analysis of specific form factors, but the underlying physics remained elusive. 

One should not expect heavy quark symmetry to apply to D meson decays, since no quark 
in the final state is heavy enough compared to the QCD scale. Nonetheless, if regarded as 
a constituent of hadrons, the strange quark has an effective mass of about 1/2 GeV, which 
is large enough that one might at least expect some vestiges of heavy quark symmetry to 
apply. In the present paper we take this point of view. We compare the predictions of heavy 
quark symmetry for the decays D — > Klv and D — > K*lu with experiment, extracting the 
nonleading form factors which account for the corrections. From these we attempt to isolate 
the physics responsible for the violation of the symmetry. The question of whether heavy- 
quark symmetry applies in this case, and even in the case of Ki 3 decays, was raised some 



time ago [|H],[T7|]. Others have shown that the leading-order predictions of heavy-quark 
symmetry are not sufficient to explain the branching ratios in semileptonic D decays. 

In brief, we have found that both the perturbative QCD and 0(l/m s ) effects are impor- 
tant. 

There is an overall suppression of decay rates to the Klv and K*lv final states due to 
perturbative QCD corrections |],§,|T{|[2(J|21|. The 0(l/m s ) corrections add together such 
that they tend to cancel this suppression in the case of Klv decay. In the case of K*lv decay, 
however, the 0(l/m s ) corrections tend to cancel each other, leaving the net suppression from 
the perturbative QCD correction. 

Since the total semileptonic decay rate for D mesons is close to that expected for free 
charmed quarks, there must be important final states besides Klv and K*lv. (We recognize 
that a recent search for such final states has yielded only upper limits P2"| , and that one 
reasonably successful description of the inclusive semileptonic rates in terms of known final 
states has appeared f23|.) We are then led to conclude that if there are indeed "missing" 
final states, they are produced primarily at the expense of the K*lv channel, which would 
help us to identify their nature. 

A previous analysis of D semileptonic decays in the same spirit as ours has appeared 



24| . Our analysis is somewhat more explicit about the sources of heavy quark symmetry 
breaking, and differs in conclusions. 

We have organized this paper as follows. We begin by describing the predictions of heavy 
quark symmetry in Section ^. After a discussion of lifetimes and branching ratios in the 
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Table 1: Kinematic factors in heavy meson semileptonic decays. 



Decay 




f max 


B -> Dev 


0.35 


1.59 


B -> D*eu 


0.38 


1.50 


D — > Ke^ 


0.27 


2.02 


D -> Jfei/ 


0.48 


1.28 



data and the leading-order heavy quark theory (Section ||), we use B decays to determine 
the Isgur-Wise function £ and the D — > K*lv differential branching ratio to determine the 
size of the perturbative QCD corrections in D decays (Section |J). We then apply the theory 
to D decay form factors in Section |5|. The implications of this analysis for B decays is 
discussed in Section |. Section [7| concludes. 



2 Heavy Quark Symmetry 
2.1 Leading-order Results 

In this subsection we discuss both the formalism and leading-order results of heavy quark 
symmetry. In the following subsections we will discuss both perturbative QCD [0(a s )] and 
finite-mass [0(l/m q )] corrections. 

It is most convenient to talk about heavy quark processes in terms of velocities instead of 
momenta. The most general velocity-dependent form factors for the decay of a pseudoscalar 



meson H to a pseudoscalar (vector) meson h (h*) |25| are: 



(h(v')\Vp\H(v)) = v^[(+W(« + 4 + ^H(«-4]. 

(h*(v',e)\V^\H(v)) = Wm H m h *iv{w)e^ape* v v' a vP, 
(h*{v',e)\A^\H{v)) = y/m H m h . [UH(w + l)e% 

-U 2 ( w ) e * ■ vv n - U 3 (w)e* ■ wjj , (1) 

where e is the /i*'s polarization vector and w = v ■ v' — {m H + rn 2 h(t) — g 2 )/(2m^m ft (,)). 
Throughout this work we denote the mass of the heavy meson H (h) by m# (m^) and the 
mass of the corresponding heavy quark by rriq (m q ). The allowed kinematic range is given 
by 

1 < w < w max , (2) 

where 

«W = ~ (z& + (3) 



2 V zW. 

and z^*> = m h (*) /mjj. Table [l] gives numerical values of z^*> and w m3iX for B and D semilep- 
tonic decays. 
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In the absence of finite-mass and radiative corrections 



(4) 



and 

£- = U = o, (5) 

where is the Isgur-Wise function. This symmetry is broken by both radiative and 

finite-mass corrections. 

The Isgur-Wise function satisfies the zero-recoil condition 

£(1) = I- (6) 
It has become standard in the literature to characterize its behavior for w near unity by 

£(«;) = (7) 

Common parameterizations of the Isgur-Wise function include the monopole, dipole, or more 
generally, the n-pole function 

£H 



(i + fOo-i) 

Another possible form is an exponential 

= exp —p 2 (w — 1 



(9) 



Other common parameterizations [25] fall in between the exponential and monopole 
forms. Since Eq. ([]) is the n —>■ oo limit of Eq. (H), we use the monopole and exponential 
forms to represent a reasonable range of forms for the Isgur-Wise function. 



2.2 Perturbative QCD Corrections 

The perturbative QCD [i. e., 0(a s )] corrections to the heavy quark limit are calculable and, 
as such, have been extensively studied [|ll[^, [T9| . ^0| .|2ll. There are two approaches to calculating 
the perturbative QCD corrections in the effective theory. The two schemes involve different 
schemes for matching the effective theory on to the full theory. One approach is to assume 
the scales m q and mq are well separated and then do matching both at m q and nig A 
second approach is to do the matching at some intermediate scale \x |19j . Hybrid approaches 
have also been developed [p(i|,|2Tl . 



For B — ► D decays these corrections are well defined. For D — ► K decays, however, the 
situation is somewhat ambiguous. Perturbation theory makes sense at the scale m c , but, 
whichever approach we use, we inevitably refer to the quantity a s (m s ). There are several 
ambiguities in using a s (m s ~ 500 MeV): There is a strong dependence on the precise value 
of m s . Even given a specific value of m s , there is a large experimental uncertainty in the 
precise value of a s (m s ). Finally, the value of a s (m s ) is large, so higher-order corrections can 
be substantial. We will have to be satisfied with finding the rough effects of perturbative 
QCD corrections for D — > K decays. 
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Table 2: Perturbative QCD correction common to all form factors in D — > ifW decays. 





^QCD 


1.00 


1.269 


1.20 


1.227 


1.40 


1.189 


1.60 


1.155 


1.80 


1.124 


2.00 


1.096 



Table 3: Perturbative QCD functions for D — > K decays 



w 




P- 


1.00 


-0.982 


-0.260 


1.20 


-1.068 


-0.254 


1.40 


-1.153 


-0.248 


1.60 


-1.238 


-0.241 


1.80 


-1.320 


-0.235 


2.00 


-1.401 


-0.229 



In the spirit of looking for the rough effects of the perturbative QCD corrections, we 
choose to use the simplest set of radiative corrections that contain the basic physics of the 
process. Therefore, for this work we use the radiative corrections given in Ref. [pO| . This 
work is a hybrid method where the matching is done at an intermediate scale and the leading 
logarithms are summed. The result is 



where 



Xqcd = Z IR {w) 



H Piiw) 



71 



a 



-6/(33-2n/) 



a s (m n 



(10) 



and Cj = {1, 0, 1, 1, 0, 1} for i = {+, — , V, Ai, A 2 , A 3 }. The mass scale \i is chosen at some 
intermediate scale between m q and tuq. See Ref. |20| for the definitions of the functions fa 
and Zir. Tables |2|-f| give numerical values of the relevant functions. 

We are now left with the ambiguity in choosing /j,, the scale at which to do the matching 
of the effective theory to the full theory. The depth of this problem is illustrated in Fig. [I], 
where we have plotted a s over the range in question. For the B — ► D case Neubert |2(]L|21[ 
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Table 4: Perturbative QCD functions for D — > K* decays 



w 


8v 


3 a 


3a 
ma 2 


3 a 


1.00 


-0.315 


-1.648 


-1.198 


-1.229 


1.05 


-0.333 


-1.644 


-1.180 


-1.231 


1.10 


-0.350 


-1.641 


-1.162 


-1.233 


1.15 


-0.368 


-1.638 


-1.145 


-1.237 


1.20 


-0.386 


-1.637 


-1.129 


-1.240 


1.25 


-0.404 


-1.636 


-1.113 


-1.245 


1.30 


-0.422 


-1.636 


-1.098 


-1.249 



has used the summation of the leading logs as a guide for picking the appropriate scale. To 
do this, notice that, under the renormalization group p|,|3|,[T9|j, 



1 + c^Q) ln QL a {m Q ) _^ ( a s (m Q 
7T a s (m q ) \a s (m, 




6/(33-2?!/) 

(12) 



We can then use this to choose the scale \i such that the two expressions are equal. Following 
this prescription for D — > K decays gives a s (fi) = 0.73. The uncertainty in a s {m c ) induces 
an uncertainty of 0.13 in this quantity. The uncertainty due to the other effects mentioned 
above is hard to estimate, but it is certainly sizable. Later we will use the data to estimate 
the value of a s (fi) appropriate forD^K decays. Note that the relevant term in Eq. ( |10|) 
is a s (/i)/7r, so the corrections due to an a s (/x) of 0.73 are 0.73/3.14 ~ 25%. 

2.3 Finite-mass Corrections 

The leading-order heavy quark results given in Eqs. (H) and (|5]) can be easily obtained using 
the trace formalism ||27|| . In this formalism 



(h^(v')\T\H(v)) = -^m hW m H ^{w)Tr[h^TH}, (13) 

where 

1 + i 

h = -^75, (14) 

h* = -p^, (15) 

h = 7o h^ (16) 

and e is the vector meson's polarization vector. H is defined similarly. In meson decays 
the current V can be either or A^. This formalism is the most convenient for considering 
the leading-order finite mass corrections There are subleading terms of 0(A/m q ) 
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and (9(A/mg), where Q (q) is the heavy quark associated with with the meson H (h). To 
0(l/m Q(q) ), ' 

A = m-jj — rriQ = m-g — m q . (17) 

The subscripts H and are barred to indicated that there is no distinction between the 
pseudoscalar and vector masses to this order. We take this average mass to be the spin 
average, i. e., = (3rrih* + m^/A. The quantity A is a fundamental parameter of QCD 
which is not calculable in perturbation theory. The experimental determination of this 
quantity is very important. Quark models suggest A = 300 MeV. With this choice of A, we 
have nib = (3m#* + m£)/4 — A = 5.02 GeV, m c = (3mo» + m£>)/4 — A = 1.67 GeV and 
m s = (3rriK* +m^-)/4 — A = 0.49 GeV. (A value A = 0.5 GeV, obtained from QCD sum rules 
28|, would imply that the quark masses are lighter, and an expansion in powers of A/2m s 



would not be valid.) In general there are corrections to infinite-mass limit of both order 
1/mg and l/m q , but we keep only the dominant (i. e., l/m q ) terms. In this approximation 
Eq. (|13"D becomes M 



(h { *\v')\T\H(v)) = y/m hW m H ^(w)Tr[h^TH] 
. A 
2m„ 



+ 



+ <) + £MK - <)J Tr [hV-frH 
- [if)+(w)(l + w)- CH(1 - w)] Tr [^h&Y rH ] } • 

To this order there are four new undetermined functions ^1,2,3,+ - These dimensionless func- 
tions are related to the dimensionful functions defined in Ref. [[| by 

Xl,2,3 = ^1,2,3 (19) 

and _ 

^ uke = (20) 

The superscript "Luke" is to distinguish the £ + used in Ref. ||, which is one of the four 
unknown functions introduced at this order, from our £ + , defined in Eq. ([I]). The functions 
ipi and ^3 satisfy the constraint || 

M™ = !) = = 1) = 0. (21) 

This constraint, which has come to be known as Luke's theorem, has been shown to follow 
from the Ademollo-Gatto theorem P9"| . 

Eqs. fll3|) and (0) can be understood schematically via Figs. ^| and |3|. In Fig. || the decay 
of the meson H is described by the decay of the heavy quark Q via the current V. This is 
represented by the trace in Eq. fll~3|). The light degrees of freedom (the dashed line) factor 
out of the trace to become the Isgur-Wise function £. At sub-leading order the current V is 
modified and the heavy quark can interact with the light degrees of freedom. The interaction 
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of the heavy quark q with the light degrees of freedom is represented schematically by Fig. [3[ 
The heavy quark q can interact in a spin-indepent way yielding the same trace, but a new 
undetermined function (ipi) for the light degrees of freedom. There can be two different 
spin-dependent interactions. The first has the spin of the heavy quark (a At!/ ) interacting with 
the light degrees of freedom and the velocity of the parent quark, yielding the function ip2- 
The factor (1 + f)j1 represents the propagation of the heavy quark. Note that 

i±£h«=/iW (22) 

The second spin-dependent interaction corresponds to the hyperfine interaction in the quark 
model. It is represented by ip^. The modification of the current V produces two new trace 
structures, but can be shown || to only require the combination of the original Isgur-Wise 
function and one new unknown function 

The net effect on the form factors is as follows [j30fl : 

£+ = £ + ^ [A - 20 - l)Va + 6^ 3 ] (23) 
£- = ^-[(2-^+0 + 1)^+] (24) 
iv = £ + 2^-K + ^- 2 ^ (25) 



2m q 



"' ^ + ^1-2^3 



w + 1 



(26) 



U = ^-[-^ + 2^2 + ^+] (27) 

u = e + [^1 - 2^ 2 - 2^3 + v+] (28) 

2.4 Momentum-dependent Form Factors 

Experiments measure the momentum-dependent form factors defined by [cf. the velocity 
dependent form factors defined in Eq. ([!])] 

(hip'W^Hip)) = f + (q 2 )(p + p')» + f-(<l 2 )(p-p')» 
(h*ip',e)\V^\H{p)) = 



{h*(p',e)\A^\H(p)) = (rriH + m^A^y^ 

J2\ 



■U<f, , , .. 

'e -p)(p + p) M - 



,e*-p)(p-pV (29) 
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The form factors /_ and A 3 lead to contributions proportional to the electron mass 
and are thus unmeasurable. The four measurable form factors are related to the velocity- 
dependent form factors defined in Eq. ([I]) as follows: 



z - 1 

7+1 



M = - — ; {w + l)£ Al , 



A; 



m h + m H 
m h + m H # 



2y/m h m H 



[z*U 2 + U s ]- 



and 



m h + m H 



(30) 

(31) 
(32) 

(33) 



2.5 Heavy Quark Symmetry to G(a s ) and 0(l/m q ) 

Considering both subleading and perturbative QCD correction terms, Eq. (j30l) becomes 



z + 1 
2^ 



0/-X"qCd£ + 



A 

2m„ 



z - 1 

ITT 



where 



and 



a 



i + — A 



7T 



2+1 
2 - 1 



* ; = Vi - 2(w - 1)^2 + 6^3 + " — t(™ + 1)^+- 

z+1 

Proceeding similarly with the other form factors, we obtain 



M = ■ [w + 1) 

m h * + m H 



where 



and 



1 + 

7T 



m h * + m H 
2y/m h *m H 



* Al = - 2^3 ; 

A //w-1 



gu 2 ^qcd£ + 



where 



and 



2m 9 VV w + 1 
l + ^( z *P A2+ P A \ 

7T v ' 
U% EE Vl + 2(^* - 1)V> 2 - 2^3 + (2* + 1)^ 



- 2* k + 



(34) 

(35) 
(36) 

(37) 

(38) 

(39) 
(40) 

(41) 
(42) 
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Table 5: Lifetimes, branching ratios (Ref. pi]), and decay rates of D mesons. 



Quantity D° D + 

Lifetime (10~ 13 s) 4.20 ± 0.08 10.66 ± 0.23 

B(D — > Xe + v e ) (%) 7.7 ±1.2 17.2 ± 2.9 

T(D -»• Xe + v e ) (10 11 s" 1 ) 1.83 ±0.29 1.61 ± 0.27 

T(D ->• [K + K*]e+u e ) (10 11 s- 1 ) 1.33 ± 0.20 0.99 ±0.12 

[K + K*)/X ratio 0.73 ±0.15 0.62 ±0.13 



and, finally, 



where 



and 



V 



m h * + m H 



a v X QCD ^ ± 7— (e ± * 



2m n 



a v 



1 + 

7T 



3 • 



(43) 



(44) 



(45) 



The functions ^//,Ai,a 2 .v represent the contribution to the form factors due to the unknown 
functions ipi. Although there are four measured quantities and four unknown functions, there 
is still a relation among the form factors: 



(46) 



3 Lifetimes and branching ratios 



In this section we shall discuss the total D meson lifetimes, the inclusive semileptonic decay 
rates, and the rates for decays to exclusive final states. We shall also compare exclusive 
decay rates with results of specific models @, 10, 11, 12, 13, 14, 15], to show that many such 
models encounter problems in describing the data. 



3.1 A Brief Tour Through the Data 

The charged and neutral D mesons have different lifetimes and semileptonic branching ratios 



31 1 , as shown in Table [5|. However, the semileptonic decay rates are compatible with one 
another. Their average is T(D — > Xe + v e ) = (1.71 ± 0.20) x 10 11 s _1 . We shall return to the 
last two lines of Table | presently. 

The experimental data on specific semileptonic decay channels are summarized in Table || 
The data we use have been summarized in Ref. ||32|| . Original sources for the branching ratios 



are Refs. ||,||,||,||,||^|3|,|9|,|g,|T],|4|,|3 
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Table 6: Semileptonic branching ratios and decay rates of D mesons to Ke + v e (Ref. 



and K*e + u, 



Decay 



Ref. 



Branching ratio (%) Rate (10 10 s 1 



D c 



K~e + v e Mark III |33 



E691 34 



CLEO §5 



ARGUS [0 



E653 |7 
Average 



3.4 ±0.5 ±0.4 
3.8 ±0.5 ±0.6 

3.8 ±0.3 ±0.6 

3.9 ±0.2 ±0.7 

2.5 ±0.4 ±0.5 
3.3 ±0.4 a ) 



7.9 ±1.0 



D+ -> fsT°e + z/ e Mark III |3| 



E691 H 
Average 



6.5 ±1.6 ±0.7 
6.1 ±0.9 ± 1.6 
6.3 ± 1.3 



5.8 ± 1.2 



Average 



7.0 ±0.8 



D c 



K* 



e v P 



Mark III |23] 



CLEO m 

Average 



3.5 ± 1.5 6) 
1.7±0.8 c) 
2.2 ±0.7 



5.2 ± 1.7 



D+ _> AT*°e + z/ e Mark III [23| 



E691 H 
WA82 H 
ARGUS [|| 
E653 e > [K2| 



Average 



4.2 ± 1.6 d ) 
4.4 ±0.4 ±0.8 
5.6 ±1.6 ±0.9 
4.2 ±0.6 ± 1.0 
3.25 ±0.71 ±0.75 
4.1 ±0.5 



3.9 ±0.5 



Average 



4.0 ±0.5 



a ^We became aware of the measurement of BR(D° — > A!" = 4.0 ± 0.8 ± 0.8% p3 

after the completion of this work. This value is not included in the average. 

b ) Based on quoted B(D -> [KTt]e + v e ) = (4.4±J;g ± 0.6)% 

multiplied by resonant fraction 0.79lo;i7^o!oi- 

c ) Based on quoted ratio B(D° -> K*-e + 'v e )/B(D° -> A"~e+z/ e ) 

= 0.51 ± 0.18 ± 0.06 times our average for B(D° — > A"~e + z/ e ). 

d ) Based on quoted -> [i?7r]e+z/ e ) = (5.3tJ;? ± 0.6)% 
multiplied by resonant fraction O^i^io;^. 

e) Based on decay £>+ — > K*°n + i' l 
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Table 7: Ratios B(D -> K*e + v e )/B(D -> ir e + z/ e ). From Ref. p 



Ref. 



Ratio 



Mark III @ 
CLEO H 
ARGUS [BHI 



E691 |? 



Average 



o.8ol°;i 

0.51 ±0.18 ±0.06 
0.55 ±0.08 ±0.10 
0.55 ±0.14 
0.57 ±0.08 



The decay rates to exclusive final states for neutral and charged D mesons are consistent 
with one another. The averages over charged and neutral decays for each channel are also 
quoted. If we sum over decays to Ke + v e and K*e + i/ e , we obtain the rates shown in the fourth 
row of Table |^. These rates fall short of the total semileptonic decay rates, though not with 
overwhelming statistical significance, as shown in the last row of Table |5|. The average of the 
two numbers in the last row of Tabled is 0.67±0.10. The Cabibbo-suppressed exclusive final 
states should account for about 5% x (relative phase space) ~ 7% of the shortfall, and the 
nonresonant fraction of the decay mode D — > [Kir\e + v e should account for another 6%. 
This still leaves about (20 ± 10)% of the inclusive semileptonic decays unaccounted for. 

The ratio Y(D — ► K*e + v e ) /T(D — > Ke + v e ) is measured to be considerably smaller than 
that predicted in most models. Experimental values for this ratio are collected in Table [7[ 

The predictions of a number of models for the ratio in Table [7] are reviewed in Ref. . 
These predictions are typically 0.9 or higher [|9], pi], [11], [12] except for the more recent lattice 
calculations [13, 14|, which are still quite uncertain. One can see a residue of the heavy-quark 
limit in these predictions. In the limit of very heavy initial and final quarks, a counting of 
spin degrees of freedom would then lead to a Vlv/Plv ratio of 3, where V and P stand 
for the ground state vector and pseudoscalar mesons [ffj. This is not far from the actual 
situation in the decays B — > D*lu and B — > Dlv. 

If the recoil of the final quark could be neglected, only the P and V mesons would be 
produced in the final state, so that P would account for 1/4 and V for 3/4 of the hadrons 
in the semileptonic final state [|I|]. By comparing the average inclusive semileptonic decay 
rate T(D -> Xe + u e ) = (1.71 ± 0.20) x 10 u s -1 with the averages in Table g we see that K 
accounts for (7.0 ± 0.8)/ (17.1 ± 2.0) = 0.41 ± 0.07 of the D semileptonic decay rate, while 
K* accounts for (4.0 ± 0.5)/(17.1 ± 2.0) = 0.23 ± 0.04 of that rate. The discrepancy with 
respect to the heavy-quark limit thus is much more marked for the K* than for the K. 

Our task is then to account for departures from the heavy-quark limit which lead to a 
modest enhancement of the rate for D — > Klv and a substantial suppression of that for D — > 
K*lv. We have mentioned in the Introduction that a combination of two effects, perturbative 
QCD and 0(l/m s ), seems to be responsible. Our subsequent discussion attempts to put 
this claim on a quantitative footing. 

Additional discrepancies with respect to the heavy-quark limit in charmed meson decays 
show up in the behavior of individual form factors [119. These are discussed in Section [5]. 
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3.2 Comparison with Leading-order Heavy Quark Theory 



It is convenient to have expressions for the widths in terms of the velocity-dependent form 
factors. Defining 



G 2 F \V Qq \ 2 m%z(*y 



48vr 3 

we find the following expressions for the branching ratios: 



dT(H -> hlv) 
dw 



T (z + l)Hw 2 -lf/ 2 



z-l 12 



(47) 



(48) 



dT T (H -> h*lv) 
dw 



2r (l + z* 2 - 2z*w)(w + l)(w 2 - I) 1 ' 2 x 

w-i)e v (w)+( W +i)e Al (w) . 



(49) 



dw 



r (w + 1)> 2 - 1) 1 ' 2 x 



(50) 



These expressions agree with those found by Neubert in Ref. ||. 

In order to illustrate the problems encountered by heavy quark symmetry for D semilep- 
tonic decays, we adapt the above relations for the decays D — > Ke + v e and D — > K*e + v e . 
We use the relations (f|) and @ of Section [24] , and predict rates as a function of as 
described in Section |2^. The results are shown in Fig. [|. 

We see that in the lowest order of the heavy-quark theory, in the absence of perturbative 
QCD corrections, the rate for the decay D — ► Ke + v e is not badly predicted, but that for 
D — > K*e + v e is about a factor of two above the experimental value. The K* / K ratio is 
predicted to lie in the range of 1.2 to 1.7 for the monopole form factor and 1.4 to 2.2 for the 
exponential form factor, in contrast to the experimental value of 0.57 ± 0.08. Perturbative 
QCD corrections lead to an overall decrease in the predicted rates and cannot account for 
the K* /K ratio. As we shall see, spin-dependent 0(l/m s ) corrections are needed for that 
purpose. 



4 Extracting Parameters from the Data 

4.1 Determining the Isgur-Wise Function from B Decays. 



We use both total and differential branching ratios in B decay to determine £(w). We have 
used the expressions ©-(W) with H = B°; h, h* = D+, D* + ; and |Vb„| = \V, h \ = 0.041 " 
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'Qq\ — I Vcb 

with the leading-order plus perturbative QCD form factors in Eq. fllOD in comparison with 
the total branching ratios given in Ref. to obtain the results shown in Table |8|. 

We can also use the B — ► D*lv differential branching ratio to determine p. This process 
has been considered in detail by Neubert M. Our fit to the ARGUS data is shown in 
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Table 8: p values obtained from various fits. 



Fitted value of p 
Decay Monopole Exponential 

Integrated Branching Ratios 
B — ► Dev 1.02±g;f| 0.93 ±0.16 
B^D*eu 0.97 ±0.21 0.91 ±0.17 

Differential Branching Ratio 
B^D*eu 1.11 ±0.50 1.02±g;^ 
Average 
1.00 ±0.15 0.93 ±0.10 



Fig. ^. As can be seen in Table our determination of p is dominated by the integrated 
branching ratios. The plot of the differential branching ratio serves as a visual check of 
our results. We have chosen the ARGUS data as representative; other data exist JRJ. The 
Isgur-Wise function can be easily extracted from the data by plotting the square root of the 
differential branching ratio divided by the the factors that multiply the Isgur-Wise function 
in the symmetry limit. This procedure can be improved by also dividing out the perturbative 
QCD corrections to the £a x piece, which dominates the rate. 
We then plot |V c6 | [(1/G) dr/dw] 1/2 , where 

G% K ^ q / — -, o r 4u> ^ Jf2 \ . ! -i *\2 



G = \V Q f^m 5 H z* 3 V^l(w + l) 2 X 2 QCD a 2 Al [^-(1 - 2wz* + z*<) + (1 - z* 

(51) 

which amounts to taking the Wq >9 — > oo and a s — > limit of Eqs. ( |4*9"D and (|5H| ) and then 
subtracting the perturbative QCD correction to the dominant part of the rate. At this level 
of approximation, [(1/G) dr/du;] = £,(w) 2 . Furthermore, the zero recoil condition 




(52) 



w=l 



holds including corrections of 0(a s ) and 0(l/m q ). 

In Refs. [jBlpEijf the quantity [(1/G) dT/dw} 1 ^ 2 was labeled simply £(w), but we prefer the 



present notation, as it makes the assumptions involved explicit. Instead of trying to fit for 
\Vcb\, we have used a fixed value (0.041) because we are only interested in the overall shape 
of the form factor. In Ref. ||, values of p = 1.14 ± 0.23 and p = 1.07 ± 0.22 were obtained 
using the pole form £(u>) = [2/(u> + l)] 2p and the exponential form Eq. (|^), respectively. A 
monopole form (Eq. (H) with n = 1) was used in Ref. [0 to obtain p = 1.26 ± 0.19. 

Both methods show a systematic dependence of the fitted value for p with the functional 
form. Fig. |^ shows that the two fitted functions look very similar with the different values 
of p. The precise value of p is important for determining \V c b\ by extrapolating to w — 1. 
We are not attempting to do this here, however. The functional dependence will have some 
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impact if we want to extrapolate beyond the end of the data (w ~ 1.4). This is the case for 
D — > Klv decays, where w max = 2.0. To minimize the systematic error from the functional 
forms, we will use p = 1.0 with the monopole form and p = 0.93 with the exponential form. 

Analyticity arguments suggest that the Isgur-Wise function should have a radius of con- 
vergence of approximately unity, i. e., it should have a pole near w = —1. This corresponds 
to a threshold for particle production near q 2 = (mi, + m c ) 2 . The n-pole function has a pole 
at w — 1 — n/p 2 . For n = 1 and p ~ 1 the pole is near w — 0. The monopole description 
can only be valid if there are a series of poles at w < — 1 which somehow conspire to mimic 
a monopole description with a pole near zero. This seems unlikely, so we conclude the de- 
pendence on the functional form of the Isgur-Wise function is even less than the difference 
between the monopole and exponential forms. In the remaining plots we use the exponential 
form with p = 0.93. 



4.2 Determining a s (fi) from D decays 



The condition in Eq. ( p2|) has been used to determine the unknown parameter \V c b\ from B 
decays. Although the corresponding parameter in D decays, |V^ S |, is well known a priori, we 
have argued in Section that the appropriate value of a s (p) is not. In this section we use 
the D — > K*lu differential width extrapolated to the point w = 1 to determine a s (p). Not 
only do l/m s corrections vanish at this point, but the QCD corrections to the dominant 
form factor are the largest of any of the QCD corrections, as can be seen from comparing 
the values of with the other /^'s in Tables || and f|. 



The published data j|2], (48], f4£j includes only raw data, i. e., including detector efficien- 
cies, compared to Monte Carlo simulations used to determine the form factors. We have 
turned this around to determine a net efficiency at each point and then determine the net 
distributions. The experimentalists themselves could do a much better job of this, but our 
method will have to suffice until they publish their own differential branching ratios. 

In Fig. ^ we have plotted the average of the two available data sets along with our fit 
used to extrapolate to w — 1. The fitted function is an exponential. The data is very flat, 
so the dependence of the extrapolation on the functional form is slight. The fit yields 



' i dr\ 

aAl ~G ~dw j 



0.457 ±0.024, (53) 



w=l 



which gives a s (p) = 1.03 ±0.11. This value is considerably larger than the guess we had in 
Section fT2| > but the uncertainties in that guess are very large. To get a consistent picture 
of D decays using only first order symmetry breaking corrections we have to accept a large 
value of a s (/i). The corrections are then 0(a s (p)/7r) ~ 33%. 



5 Form Factors 

5.1 Comparison with Experiment 



All four D — > K form factors have been measured. Three groups [[!3|,|33|,|3l|] have measured 
the q 2 dependence of /+. They find that the data are consistent with the monopole vector 
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Table 9: Form factor measurements for D — > K*lv. 



Experiment 


A 2 (0)/A 1 (0) 


nO)Mi(O) 


AM 


A 2 (0) 


7(0) 


E691 


0.0 ±0.5 ±0.2 


2.0 ±0.6 ±0.3 


0.46 ±0.07 


0.0 ±0.2 


0.9 ±0.3 


E653 


0.82 ±0.23 


2.0 ±0.3 ±0.2 


0.49 ±0.07 


0.40 ±0.14 


0.99 ±0.27 


Average 






0.48 ±0.05 


0.27 ±0.11 


0.95 ±0.20 



dominance form 

/+(? 2 ) = i /+ il 2 - (54) 
1 — q z /mj 

We use the world average value given by Stone 

ro, = 2.0lg|. (55) 

The normalization constant can be obtained from integrating the branching ratio. Here, too, 
we use the average value cited by Stone |32|: 



/ + (0) = 0.69 ±0.04. (56) 

The three measurable form factors in D — > iCZv semileptonic decay have also been measured 
TJ.[I5|.[ID|]. Unfortunately, these measurements assume the g 2 dependences follow single-pole 



vector dominance so that 

Ai, 2 (0) 
1 — q 2 jm\ 

and 

V(0) 



= rl^ (57) 



^ ) = 1 (5 8 ) 

1 — q z /m v 

where ttia and my are assumed to be the lowest mass cs states with J p = 1 + and 1~, i. e., 
ttla = 2.5 GeV and my = 2.1 GeV. The experimentalists then fit for the ratios A2(0)/Al(0) 
and V(0)/Ai(0). ^(0), A\(0) and V(0) then may be determined from the total branching 
ratio, which is dominated by Ai(0). 

Table |9| displays the values obtained by the two experiments. The measurements of V(0) 
and Ai(0) are quite consistent with one another. The values of A 2 (0) differ from one another 
by slightly less than two standard deviations. 

It should be pointed out that there is some confusion stemming from the way these values 
are reported. The numbers V(0), ^4i(0) and A 2 (0) have been interpreted as measurements 
of the these form factors at a single point (q 2 = 0). In fact, both experiments do their fits 
to Eq. (|5"7| ) and Eq. (|S8) over the entire range of q 2 . V(0), Ai(0) and A 2 (0) are merely the 
constants in Eq. (|57|) and Eq. (|58"D . 

In Fig. we compare the measured form factors with our predictions from heavy quark 
symmetry. We see that the deviations from the leading-order heavy quark predictions are 
about as large as was expected: ~ A/2m s ~ 30%. The rates depend on the square of the 
form factors, however, so a deviation in the form factor of 25% can lead to a difference in rates 
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by a factor of 2! This helps explain the discrepancies seen in Section |3|. Note that, although 
the D — > Klv rate was closest to the leading-order prediction of heavy quark symmetry, it 
has the largest deviation in the form factors. The l/m s effects cancel the a s effects for this 
decay, while they add coherently in D — ► K*lv decay. 

5.2 Determining A from V(w = 1) 

It follows from Eqs. ©, (|l]) and © that 



V(w = 1) = ^j=== (ayXqcu + ■ (59) 

A measurement of V(w = 1) is the theoretically cleanest way to measure the fundamental 
quantity A. Even in D decays, where the size of is in question, this relation is useful 

because the dependence on a s (n) is relatively weak. Using the measured V from the previous 
section, this leads to 

— — = -0.05 ± 0.239, (60) 
2m s y ' 

and 

A = -44 ± 297 MeV. (61) 

This error is clearly too large to be truly useful. We have the further caveat that the 
extrapolation of V(w) to w — 1 will only be reliable once the w-dependence of the form 
factors is measured. We hope this quantity will be measured with smaller uncertainties in 
the future. 

The above measurement of A is subject to the various uncertainties involved in applying 
heavy quark symmetry to the strange quark. It is much more important to measure A/ (2m c ), 
where heavy quark symmetry is on very sound theoretical ground. CLEO 5(| has recently 



published a first measurement of V(w = 1) in B — ► D*lu decays. They obtain two different 
values from two different fits: 

K(m=1)H °->™9±0.12(fita) 
1 1.19 ±0.57 ±0.15 (fit b) 

which we average to get V{w = 1) = 1.05 ± 0.5. Using mu = 5.28 GeV, rrih* = 2.01 GeV, 
Xqcb = 1-11 and ay = 0.97 in Eq. (|59D leads to 

A -0.14 ±0.45. (63) 



2m 



The uncertainties are so large that it doesn't make sense to extract A. It is important that 
this quantity become better measured in the future. 
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5.3 Extracting the Subleading Effects From the Data 



In Fig. H we have assumed that the differences between the predictions and the experiments 
shown in Fig. [I] are completely due to subleading terms. We have then subtracted out the £ 
dependent effects to obtain \l/ f t A u A 2 ,v ■ This procedure depends on the value of A, so we have 
included results for both A = 200 and 400 MeV. The size of the experimental uncertainties 
depends on A because we have scaled out A/(2m s ). 

Several things are apparent from Fig. || First, the dimensionless functions all come out 
to be 0(1). We also see that the prediction ^1% = tyy is better satisfied for the smaller value 



of A. This is consistent with the results in Section |5.2| . We also see that the corrections to 
the pseudoscalar channel embodied by ^ are substantial and positive. The net effect of the 
corrections to the vector decays is small and negative. 

For concreteness, we can attempt to estimate the individual functions ipi under certain 
assumptions. First, we choose the smaller value of A in Fig. |8|. Using the definition of ^ f 
[Eq. fl36f )1 and ipi(l) = ^3(1) = (Luke's theorem), we obtain 

= -1-15^+(1) (64) 

Comparing with the light band in the plot of in Fig. |8|, we obtain t/>+(1) ~ —1.5. Similarly, 

we use 

#a 2 (1) = -1.04^ 2 + 1.48^+ (65) 

to obtain ^(l) ~ — 2. To get information about ipi and ^3 we need to go to another 
kinematical point. It is best to use the minimum recoil point of the pseudoscalar decay, since 
it has the larger kinematic range and the kinematic dependence of the the pseudoscalar form 
factor (/) has been measured. Using Eq. (|36| ) once again, we have 

^(2) = Vi(2) - 2^(2) + 6^(2) - 1.7^+(2). (66) 

Since ^ Ai is consistent with zero, we set i[>i(w) = 2ip^{w). If we also assume i/j+(w) and 
■02 (w) are roughly constant in w, we can then obtain ^3(2) pa —0.5 and ipi{2) ~ —1. These 
estimates are very rough, owing to both experimental uncertainties and our assumptions. 



6 Implications for B Decays 

The overall QCD decrease of semileptonic decay rates into pseudoscalar and vector mesons 
found for D meson decays implies that for B meson decays similar, but slightly weaker, 
effects are to be expected. In view of the ambiguity reflected in Fig. we can make only 
qualitative statements. 

While the 0{l/m q ) effects discussed above are expected to reduce the ratio T(B — ► 
D*lv)/T(B — > Dip) from its infinite-mass limit [|IJ of 3, the lowest-order prediction of this 
ratio exceeds 3 slightly for reasonable slopes of the Isgur-Wise function [|4"7|j . The experi- 
mental branching ratios @ are B(B -> D*lu) = (4.4 ± 0.3 ± 0.6)% and B(B -> Dlu) = 
(1.6 ± 0.3 ± 0.2)%; their ratio is 2.75 ± 0.75. The error is too large to reflect the presence of 
any of the effects mentioned here. 

The sum of the semileptonic B branching ratios to D and D*, (6.0 ±0.8)%, is only 1/2 to 



2/3 of the inclusive semileptonic branching ratio [J4|, B(B — >• Xlv) ~ 11%. The remaining 
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1/2 to 1/3 of the semileptonic decays are expected to involve excited D mesons or extra 
pions. A non-trivial Isgur-Wise function (one with p > 0) is sufficient to reflect the presence 
of such excitations Jl7|]; one does not require 0(1/ m) effects to understand them. 

Perhaps the most pressing question about heavy quark symmetry and B decays is how 
much we can learn about V c b- Our analysis suggests that measuring V c b by extrapolating the 
B — > Dlv spectrum to the normalization point would be unreliable because of the substantial 
A/m c effect at that point. Fortunately, extrapolating the B — > D*lv spectrum, as was done 
in Ref. ||, avoids all complications due to A/m c effects. The question in this case, then, is 
the size of the (A/m c ) 2 corrections. We have shown in Section [)72] how the quantity A/ (2m c ) 
can be extracted from the data. Although the error we obtained is very large, the results 
favor a small value of A, which would be good news for an accurate determination of V c b if 
more accurate measurements confirm that A is indeed small. 



7 Conclusions 

We have applied a broken version of heavy-quark symmetry to the semileptonic decays 
D — > Klv and D — > K*lv, in order to study the deviations from the symmetry in an 
environment where they are particularly prominent. The price that we pay is the possible 
loss of validity of the symmetry. We have to regard the strange quark as heavy, which we 
can do if we regard it as having a constituent- quark mass of order 1/2 GeV. 

We have been able to identify several physical sources of the apparent deviations from 
heavy-quark symmetry (and of the shortcomings of early models for D semileptonic decays). 

First, perturbative QCD effects lead to an overall reduction in predicted D — > Klv and D ■ 
K*lv rates with respect to the lowest-order theory. One expects the remainder of the semilep- 
tonic decays to show up in excited kaons or K + mr final states. No evidence for these final 
states exists at present. 

Second, 0(1/ m s ) effects give different contributions to the pseudoscalar and vector final 
states. These effects add to cancel the perturbative QCD effects in the pseudoscalar final 
state. In the vector final state, however, these effects cancel one another, leaving the overall 
supression due to the perturbative QCD effects untouched. This explains the suppression of 
the vector final state relative to the pseudoscalar final state. 

We find that it is instructive to view the strange quark, under some circumstances, as 
"heavy," which allows us to gain a qualitative understanding of effects in D decays which 
are still too small to show up in the decays B — > Dlv and B — > D*lv. A study of these 
decays with increased precision would provide an excellent check of our conclusions. 
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Figure Captions 



Figure 1: Variation of a s over the region of interest for D decays based on [26 
a s (m c ) = 0.276 ± 0.014 The solid line represents the central value, while the dashed 
lines represent the one standard deviation uncertainties. 



Figure 2: Diagram for heavy meson decays at leading order. 



Figure 3: Sub leading diagram where the heavy quark interacts with the light degrees of 
freedom. 



Figure 4: Predictions of heavy-quark symmetry for decay rates of D — > Ke + v e (a,c) and 
D — > K*e + u e (b,d). Predictions without QCD corrections (i. e., a s (p) = and Xqcd = 1) 
are shown in (a,b) as functions of p for monopole (solid lines) and exponential (dashed lines) 
form factors. The shaded bands correspond to the experimental values with one standard 
deviation uncertainties. The ranges of p allowed by the fits to B decays in Section 3~T 



are shown for monopole (circles) and exponential (diamonds) form factors. Predictions for 
central values of p are shown as functions of a s (p) (holding Xqcd fixed) in (c,d), where the 
bars indicate the range of a s (p) obtained in Section |Ol 



Figure 5: Data taken from Ref. f45fl . The solid and dashed lines are fits to the exponential 
and monopole forms, respectively. 



Figure 6: Differential D — > K* distribution from Refs. [48|] and [49]. The line is a fit used 



to extrapolate to w = 1. The method of extracting this plot from the published data is 
described in the text. 



Figure 7: Comparison of the leading order heavy quark predictions with data for the four 
measured D meson semileptonic decay form factors. The solid and dotted lines correspond 
to the central value of a s (p) obtained in Section fl| and the one standard deviation un- 
certainties, respectively. The shaded bands correspond to the experimental values with one 
standard deviation uncertainties. 
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Figure 8: Values of form factors associated with subleading operators obtained from experi- 
ment. The light and dark regions correspond to taking A = 200 and 400 MeV, respectively. 
The light regions have been extended horizontally for clarity in some cases. 
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